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Abstract. We have performed a new and homogeneous analysis of all the Li data available in the literature for 
main sequence stars (spectral-types from late F to K) in open clusters. In the present paper we focus on a detailed 
investigation of MS Li depletion and its time scales for stars in the 6350-5500 K effective temperature range. For 
the first time, we were able to constrain the age at which non-standard mixing processes, driving MS Li depletion, 
appear. We have also shown that MS Li depletion is not a continuous process and cannot be simply described 
by a t~" law. We confirm that depletion becomes ineffective beyond an age of 1-2 Gyr for the majority of the 
stars, leading to a Li plateau at old ages. We compared the empirical scenario of Li as a function of age with the 
predictions of three non-standard models. We found that models including only gravity waves as main mixing 
process are not able to fit the Li vs. age pattern and thus this kind of mixing can be excluded as the predominant 
mechanism responsible for Li depletion. On the other hand, models including slow mixing induced by rotation 
and angular momentum loss, and in particular those including also diffusive processes not related to rotation, 
can explain to some extent the empirical evidence. However, none of the currently proposed models can fit the 
plateau at old ages. 

Key words. Stars: abundances - Li - Stars: Evolution - Open Clusters and Associations: Individual: 



1. Introduction 

Lithium, as well as the other light elements beryllium and 
boron, is burned at relatively low temperatures in stellar 
interiors (~ 2.5 x 10® K). As a consequence, it survives only 
in the external layers of a star and it is a very powerful 
tracer of mixing mechanisms at work in stellar structures. 

A huge amount of observational and theoretical work 
has been devoted to the understanding of Li and its evo- 
lution. In particular, during the last two decades, several 
datasets for Li in open clusters and in both Population i 
and II field stars have been collected. Focusing on Pop. i, 
observations of unevolved stars in open clusters and the 
field evidenced several features that disagree with the pre- 
dictions of the standard models. With "standard" or "clas- 
sical" we refer to those models that include convection 
only as a mixing process and do not take into account 
transport phenomena like diffusion, gravity waves, angu- 
lar momentum loss and transport, etc. The main open 
problems concerning Li evolution during the pre-main se- 
quence (PMS) and main sequence (MS) phases of solar- 
type and lower mass stars have been recently discussed by 
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Jeffries H2005f) and Randich H2005|l and are summarized 
below: 

(t) PMS: G-type stars hotter than - 5300-5400 K in 
zero-age main sequence ( ZAMS, ~ 30-100 Myr) open 
clusters (e.g., Soderblom et al. 11993 6* — hereafter S93b; 
Randich et al. 11998. Randich et al. |2001.) have Li abun- 
dances (log7i(Li)) only slightly below the initial abun- 
dance for Pop. I stars (logn(Li)o — 3.1 — 3.3, see, e.g., 
Jeffries I2005|l ; a larger amount of PMS Li depletion is in- 
stead expected from standard models. Note however that 
the predicted Li destruction strongly depends on the phys- 
ical assumptions adopted in stellar codes. In addition, K- 
type stars of presumably similar effective temperatures 
(Tcff) in clusters younger than ~ 250 Myr are character- 
ized by a wide dispersion in Li abundance whose origin is 
not yet understood, although different explanations have 
been proposed. 

(ii) MS: at variance with standard model predictions, 
the empirical evidence for clusters older than the Pleiades 
(e.g., the Hyades; Thorburn et al. —hereafter T93) 
suggests that solar-type stars do deplete their photo- 
spheric Li after having reached the ZAMS. Another puz- 
zling feature concerns the discovery of a "Li dip" for F- 
type stars around ~6700-7000 K (e.g. Balachandran ll99'51 
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and references therein): stars of the Hyades and clusters 
of similar (and older) ages within this T^s range have Li 
abundances lower by a factor 30 (or more) than those 
of stars within a Tgg range of ~ ± 300 K. Such a gap in 
the Li distribution may start appearing at ages as young 
as ~200 Myr (see Steinhauer & Delivannis 121)041 and ref- 
erences therein). 

Observations of old open clusters show that past the 
Hyades age Li depletion for solar analogs can be either 
very slow or very fast: in fact, solar- type stars in the so- 
lar age/solar metallicity M 67 are characterized by a large 
amount of dispersion in Li (see Pasquini et al. 119971 Jones 
et al. 119991 and references therein) , with about 40 % of 
the stars more Li depleted (by factors 5-10) than the re- 
maining fraction ("upper envelope"). On the other hand, 
the other old open clusters investigated so far (IC 4651 
and NGC 3680, Randich et al. 2000; NGC 188, Randich 
et al.Enm NGC 752, Sestito et al. l2004| have tight Li dis- 
tributions with abundances similar to those in the upper 
envelope of M 67 and only ^2 times lower than those of 
similar stars in the Hyades. We mention that preliminary 
analysis of the old cluster Cr 261 suggests that this cluster 
could also be characterized by a dispersion, although not 
as large as the one observed in M 67 (Pallavicini et al. 
I2005|l . A large spread similar to that present in M 67 is 
also observed for solar analogs in the field (e.g., Pasquini 
et al. I1994|l : the Sun is representative of Li-poor stars in 
the field, and its Li abundance cannot be reproduced by 
standard models. 

(in) Finally, Li depiction does not seem to significantly 
depend on the cluster metal content; both young and old 
clusters with different metallicities — varying within ^ 
±0.2 dex around the solar [Fe/H] — seem to share a very 
similar Li distribution (e.g., Jeffries & James [TI)99l Sestito 
et n]. mm^OOA ]. 

These puzzling features, not expected from standard 
models, suggest that Li depletion is not driven uniquely by 
convection, and that extra-mixing processes (or processes 
able to inhibit Li depletion during the PMS) not included 
in classical theory are at work in stars during the various 
evolutionary phases. Several efforts have been done in the 
last years aimed to understand this (these) process (es) 
and non-standard models have been developed, but so far 
the mechanism(s) driving Li depletion remain(s) poorly 
constrained. 

In the present paper we focus on point (ii) above, and 
more specifically on the detailed investigation of MS Li 
depletion and its time scales in solar- type stars. Since pro- 
posed non-standard mixing processes have different time 
scales of Li depletion, our main goal is to put empirical 
constraints on these time scales for stars with different 
temperatures/masses and thus to provide feedback to the 
models. In order to carry out such an analysis, not only 
the number of observed open clusters and of stars per clus- 
ter should be enlarged to have a good age sampling, but 
— most important — all the already available Li data for 
open clusters should be analyzed with the same method, 



to reduce the possibility of spurious results due to an in- 
homogeneous analysis. 

For this reason, we collected all the available informa- 
tion on Li in late F, G and K-type MS stars in open clus- 
ters and we performed a new analysis of Li abundances^ . 
Our final objective is to provide a homogeneous (at least 
as far as the analysis is concerned) database for open clus- 
ters, that can be used to study Li evolution and to test 
theoretical models. We stress that, whereas most of the 
papers on Li in open clusters have so far focused on the 
comparison of samples with similar ages or in a restricted 
age range, our goal here is to investigate the evolution of 
Li and its timescales from the ZAMS up to very old ages. 
The database will be also made available in the web. Note 
that the discussion in the present paper focuses on MS 
Li depletion in solar-type stars, but we included in the 
database also warmer and cooler stars (F and K-type) in 
order to provide a more complete picture for Li evolution 
in open clusters. F and K-type stars will not be discussed 
here, since they are characterized by different structures 
and evolutionary histories with respect to solar analogs 
and deserve to be studied separately. 

As a final but fundamental remark we would like to 
emphasize that, even though our method of analysis can 
be affected by systematic errors, we are interested in the 
investigation of differential depletion of Li during the MS 
and not in the absolute abundance values for each cluster: 
in this context, the importance of performing a consistent 
analysis with all the samples on the same scale appears 
even more evident. 

In Sect. 2 we present the data sample, while in Sect. 3 
the method of analysis is described; Sect. 4 includes the 
results and a discussion; some conclusions are reported in 
Sect. 5. 

2. Data samples 

Our sample includes 22 clusters (listed in Col. 1 of Table^l 
with different ages and metallicities (Cols. 3 and 8; the ref- 
erences for age and [Fe/H] are reported in Cols. 4 and 9 
respectively, while Col. 2 lists references for Li data). Cols. 
5 to 7 list the distance modulus and reddening with the 
proper reference(s). For the present analysis we considered 
stars with 4000 < Tes< 6500 K; as mentioned, the discus- 
sion in Sect. 4 will concern solar-type stars, i.e. objects 
with 5500 < Teff< 6350 K. 

The ages of the clusters range from ^ 5 Myr for 
NGC 2264, to ~ 6-8 Gyr for NGC 188, while [Fe/H] values 
are within an interval of ^ ±0.2 dex from solar. We consid- 
ered the "classical" ages for each cluster, i.e. those deter- 
mined from isochronc fitting. We would like to stress here 
that, as Table 1 shows, our assumptions for cluster ages are 
a collection of independent determinations, i.e. they were 
derived by different authors and using disparate models. 
Thus our dataset is not homogeneous as far as ages are 
concerned. In addition, we mention that for four young 

^ The database is available upon request to P. Sestito. 
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Cluster 


Rcf . 


Ago 


Rcf. 




E{B — V) 


Rcf. 


[Ic/HJ 


Ref. 


Tcff range 


No, of 


S/N 




Li data 


(Gyr) 


age 


(mag) 


(mag) 


distance, reddening 








stais ( ) 






1 


0.005 


2 


9.4 


0.071 


57 


— 0.15 


77 


4050-6350 


12(4) 


50-100 




2 
















41U0-555U 


37(8) 


25—105 


IC 2602 


3,4 


0.030 


43 


6.11 


0.04 


58 


— 0.05 


4 


4150-6450 


31(1) 


30-100 


IC 2391 


5 


0.030 


5 


6.05 


0.01 


59,60 


— 0.03 


4 


4100-6050 


6 




IC 4665 


6 


0.035 


44 


8.3 


0.18 


61 






4750-5950 


13 




NGC 2547 


7 


0.035 


45 


8.30 


0.06 


62 


sub-solar 


7 


4100-5100 


7(2) 


80-190 


a Per 


8,9 


0.050 


8 


6.06 


0.08 


63 


— 0.05 


78 


4450-6450 


39 


100-200 




10 
















4150-5400 


18 




AT o /I r" 1 

NCjC 2451 


11 


0.050 


46 


7.57 


0.05 


64,65 


— 0.01 


11 


4350-6350 


17(5) 


20-150 


Pleiades 


12 


0.07 


47 


5.36 


0.04 


66,67 


— 0.03 


78 


4100-4700 


7 






13 
















4000-6500 


95(22) 






14 
















Af\r\n /loco 

4000—4250 


4 


9u— oOU 




15 
















4050-5250 


10 


~70 




16 
















4750-5150 


9(1) 


90-110 


Blanco 1 


17 


0.1 


48 


6.9 


0.02 


68 


+0. 14 


17 


4550-5850 


17 


~30 


NCjC 2516 


18 


0.15 


49 


8.18 


0.12 


69 


~ —0.10 


18 


4810-6250 


22(3) 


10-65 


M 35 


19 


0.2 


50 


9.60 


0.225 


50 


— 0.21 


19 


4750-6250 


27 


40-160 


NGC 6475 


20 


0.25 


49 


7.35 


0.06 


70,20 


+0.14 


22 


4450-6400 


24(6) 


15-50 




21 
















5150-6200 


10 






22 
















5000-6300 


32(11) 


50-150 


M 34 


23 


0.25 


49 


8.28 


0.07 


71,72 


+0.07 


79 


4150-6450 


38(3) 


~70 


Coma Ber 


24 


0.6 


51 


4.71 


0.013 


73,64 


— 0.05 


80 


6200 


1(1) 


75-105 




25 
















4200-6300 


14(4) 


150-250 




26 
















4950-6200 


13(2) 


200 


NGC 6633 


27 


0.6 


52 


7.77 


0.017 


27 


— 0.10 


27 


4600-6500 


22(3) 


~100 


Hyades 


24,28 


0.6 


53 


3.31 


0.01 


53,64 


+0.13 


78 


4900—6200 


14(1) 


90—120 




29 
















4600-6400 


58(31) 


150—220 


Praesepc 


30 


0.6 


30 


6.39 


0.009 


64 


+0.03 


80 


4850-6450 


57(22) 


~100 


iNLrC 752 


31 


2 


54 


8.25 


0.035 


74 


+0.01 


32 


4900—0350 


5 


30—80 




32 
















5550-6300 


18(1) 


50-80 


NGC 3680 


33 


2 


54 


10.5 


0.05 


75 


— 0.17 


34 


5900-5950 


2(1) 


30-60 




34 
















5900-5950 


4(1) 


60-150 


IC 4651 


33 


2 


54 


10.1 


0.083 


76 


+0.10 


81 


5800-6250 


10(1) 


30-60 


M 67 


35 


5 


55 


9.69 


0.05 


56,39 


+0.05 


82 


5650-6200 


6 






36 
















5900-6000 


3 


20-70 




37 
















6150-6200 


3(1) 


50-100 




38 
















5900-6150 


7(2) 


~50 




39 
















5500-6150 


25(1) 


60-100 




40 
















5900-6150 


2 


100-200 


NGC 188 


41 


8 


56 


11.44 


0.09 


56 


+0.01 


42 


5750-5950 


2 


~45 




42 
















5700-6000 


11 


20-35 



* — The number in brackets in Col. 11 indicates the number of binaries present in the samples. 

References: (1) King |19ijBl; (2) Soderblom ot al. 119991 : (3) Ra ndich e t al. 119971 : (4) Randich et al. 120011 : (5) Stauffer Jlgggh (6) Martin & 
Montes (19971; (7) Jeffries et al. (2003); (8) Balachandran et al. (19881; (9) Balachandran et al. I1996|: ClOl Randich et aL^lg^; (11) Hiinsch 
et al. (20041; (12) Butler et al. (1987 1; (13) Soderblom et al. I 1993b S93b); (14) Ga rcia L opez et al. n994l : (15) Jones et al. ro96l : (16) Jeffrie s 
I1999bl : (17) Jeffries & James (1999i: (18) J effries et al. fT598ll : (19) Barrado et al. I2U01I : (20) James fe Jeffries 119971 : (2 1) Jam es et al. f2000"l; 
(22 ) Sest ito et al. (2!03l: (23) Jones et al. 119971 : (24) Sod erblom et al. (199.01; (25) Jeffries \W^^; (26) Ford et al. 120011 : (27) Jeffries ot 
al. 120021 : (2 8) Sod erblom ct al. (119'95l ; (29) T horburn et al. 119931 T 93j: (30) Soderblom et al. g9933l; (31) Hobbs & Pilach owski tl986a< : (32) 
Sestito et al.J20lH|; (33) Randich et al. J2000|; (34) Pasquini et ah 1 20011: (35) Hobbs & Pilacho wski (1986b'l; (36) Spite et al. n987l : (37) Garcia 
Lo pez et al. toSF (38) Pasq uini et al. fW^Tv . (39) Jones ot al. ('iragl; (40) Randich e t al. 120021 : (41)Hobbs & Pilach owski 119881 : (42) Rand ich et 
al. (2"003 I ; (43) St auffer et al. <ii>§7l : (44) Mermilliod et al. (1981 1; (45) Naylo r et al. JliOOl^l : (461 Hiinsc h et al . JOTS); (47) Patenaude ii07tjl : (48) 
Panagi & O'Dell ilj^^TI ; (49) M eynet et al. (ITO); (50) Sung & Bcssol (ITO); (51) Gar cia Lo pez et al. (52) S trobel OMD; (53 ) Perr yman 

et al. treSBI : (54) Friel |ii>§51 : (55) Hobbs & Th orbur n ITOlt : (56) Sarajedini et al. jii Mj ( 57) Sung et al. iggg}; (58) Whiteoak ilO^iil : (59) 
Becker & Fenkart I197q|; (6 0) Patten & Si mon 119961 : (61) Prossor (1993l|j62) Claria ifllfel : (63) Mitcheh ^601: (64) Lynga Catalogu e (198 7 
— fifth edition; Lvnga|1985l: (65) Feinstoin J^66| ; (66) Robichon ct al. (.19991: (67 1 Stauffer & Hartmann f 19871; (68) Epstein & Epstein 11985J; 
(69) Dachs & Kab us tBSST (70) Snowden ffiTO: (71 1 lanna & Schlemmer 09931 ; (72) J ones & Pressor I 19961; (73) Odenkirchen et al. (19981; 
(74) Daniel et al. ^9941 : 175) Nordstrom et al. | l9 971: (76) Anthony-Tw arog e t al. 119881 : (77) King (| lM 3i: (78) Boesgaard & Friel JT99QJ : (79) 
Schulor et al. 120031 : (801 Friel & Boesgaard 119921 : (811 Pasquini et al. 120(141 : (82) Randich et al. I2005i . 



clusters included in our sample (IC 2391, NGC 2547, a 
Per, Pleiades) new independent estimates of the ages have 
been obtained based on the position of the "Li depletion 
boundary" (LDB; see Stauffer et al. HHMI EMni Barrado 
y Navascues et al. HdM Jeffries & Oliveira EIEU) • LDB 
ages are in general higher than the classical ones: in the 
references mentioned above ages of 50 Myr for IC 2391, 
90 Myr for a Per and 125 Myr for the Pleiades have been 
reported, to be compared with those listed in Table 1. In 
the case of NGC 2547 (Jeffries & Oliveira l200^ there is 
instead a perfect agreement between the LDB age and the 
classical one ('~35 Myr). 

The classical methods for estimating the age of a clus- 
ter, i.e. the fit of the turn off or of the low MS, are affected 



by uncertainties mainly due to physical inputs adopted in 
models; on the contrary, the LDB technique is less model 
dependent. However, we adopt here the classical ages to 
be consistent with the whole sample, since LDB ages are 
available only for four clusters. 

The metallicities listed in Table were all derived 
from high- resolution data (see the quoted references), with 
the exception of those of NGC 2264 and NGC 6633 that 
were estimated from intermediate resolution spectra (King 
119931 Jeffries et al. I2002|l ; however these are the most re- 
cent and reliable estimates for the Fe content of the two 
samples. No estimate of the metallicity is available for 
IC 4665; for NGC 2547 we only indicate that the cluster 
has a sub-solar Fe content (see also Jeffries et al. I2003|l 
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since no precise estimate is present in the literature. As 
in the case of ages, we caveat that also the [Fe/H] values 
for most clusters were not derived with a homogeneous 
analysis, with exception of IC 2391, IC 2602, NGC 2451, 
NGC 6475, NGC 752, M67 for which we have derived the 
metallicity on the same scale as for the Hyades. 

Several clusters were studied by more than one au- 
thor; in these cases we identified the stars in common 
between the different samples and we selected the more 
recent observations (which have in most cases a better 
signal-to-noise ratio — S/N — and/or higher resolving 
powers — R); among the various clusters the data se- 
lection for the Pleiades deserves a slightly more accurate 
description. For this cluster we adopted the following pri- 
ority order: (1) data of Jeffries l)1999b|l : the author carried 
out observations of a sample of cool Pleiades stars in or- 
der to investigate the variability of Li i and K i and its 
possible relationship with the spread among K-type ob- 
jects. For this reason several objects were observed twice 
(November 1997 and November 1998): since the equivalent 
widths {EW) measured from the two sets of spectra are in 
good agreement (no Li variability was found over a 1 year 
time scale) we chose the value obtained from the spectrum 
with the highest S/N ratio. In one case the two spectra 
have the same S/N and thus we computed an average of 
the EWs\ (2) data of Jones et al. 1)19961 although they 
have lower S/N than those of Garcia Lopez et al. H1994<) , 
the resolution is higher); (3) data of S93b, if the spectra 
has a "a" quality (see the quoted reference), otherwise (4) 
data of Garcia Lopez et al. (19941. We preferred the "a" 
quality data of S93b with respect to those of Garcia Lopez 
et al. I|1994|l since they have in principle a slightly better 
S/N] we note however that this subsample contains only 
two stars, for which the errors in EW quoted by the two 
authors are very similar; (5) data of Butler et al. H1987|l . 
Note that S93b re-observed some stars included in the 
sample of Butler et al., but they quoted a larger error in 
EW with respect to the previous paper (and thus, proba- 
bly a worse S/N). For this reason, we retained the original 
measurements by Butler et al. I)1987|l . Finally, (6) all the 
remaining stars included in the sample of S93b. 

The number of stars and the relative range of effective 
temperatures (as computed by us, see Sect. I3.1|l for each 
sample are reported in Cols. 11 and 10 of Tabled while 
the references for the source papers are listed in Col. 2. 
Most obviously only confirmed cluster members were con- 
sidered. The number in brackets in Col. 11 represents the 
number of binaries present in the samples; double systems 
can be safely included in our analysis since in general they 
do not affect the Li distributions (but see also Sect. I4.2|l . 

We finally mention that the cluster datasets consid- 
ered here were observed during the last 20 years us- 
ing different telescopes and/or instruments; therefore, the 
various subsamples are characterized by disparate S/N 
and R. Spectral resolutions range from ~ 1 A in the 
case of NGC 2264 (Soderblom et al. up to - 

0.1 A (i? ~ 57, 000) for NGC 752 and NGC 188 (Sestito et 
al.ElSI Randich et al.EIini; S/N ratios (listed in Col. 12 



of Table nj range approximatively from ^ 20-30 — e.g., 
in the cases of Blanco 1 (Jeffries & James ^99) and of 
the old and rather distant cluster NGC 188 (Randich et 
al. .2003,1 — up to the very good quality of the spectra 
of some stars in in the Pleiades {S/N ^ 200 — 300 — 
Garcia Lopez et al. CPT), in a Per (S/N ~ 100 - 200 

— Balachandran et al. 11988,) . in NGC 6475 {S/N ~ 150 

— Sestito et al. I2003|l and other young/intermediate age 
clusters. In other words, from the point of view of the 
data quality the complete set cannot be considered ho- 
mogeneous; homogeneity would imply that all the stars 
in all the clusters should have been re-observed with a 
unique instrument and possibly in similar weather condi- 
tions. Thus, a homogeneous dataset on the point of view 
of S/N cannot be obtained, since also spectra collected 
during the same night can be characterized by discrepant 
qualities due to changes in weather conditions. 

Errors in EWs depend on the S/N ratio of the spectra, 
but also the spectral resolution affects the estimate of the 
Li EWs: as we will explain in more detail in Sect. 13.21 this 
spectral feature can be blended with a nearby Fe i line. 
When the resolution is not large enough (or in presence 
of large rotation) the two lines cannot be separated and 
we subtract the contribution of the Fe i feature to the to- 
tal EW using an analytical expression in which EW{Fe) 
is a function of S — V^. As a consequence, Li abundances 
derived from spectra with different resolving powers are 
differently affected by errors introduced with the theoret- 
ical estimate of the Fe contribution (see Sect l3.3.I)l . 

3. Li analysis 

3.1. Effective temperature 

We started our new analysis from the published EWs of 
the Li I A67G7.79A line and from BV photometry used 
in the original papers; the proper references will be re- 
ported in the web database. As already mentioned, the 
adopted reddening values towards each cluster are listed 
in Col. 6 together with the reference (Col. 7). Note that 
the E{B ~ V) values adopted for each cluster are rather 
"safe", in the sense that discrepancies with other esti- 
mates reported in the literature are very low (~ 0.03 
mag at most); this is due to the fact that the clus- 
ters considered in the current analysis are relatively close 
and well studied. In addition, none of the sample clus- 
ters is characterized by differential reddening (see the 
database by J.-C. Mermilliod and references therein — 
http : / /obs www . unige . ch / webda ) . 

We decided to use BV photometry since for the large 
majority of F, G and K-type stars there are by far more 
B ~ V observations available than for any other color in- 
dex: this allows us not to restrict the dataset too much 
(for example VI photometry is in general available only 
for cool stars). In addition, and most important, it has 
been shown (Soderblom et al. I1993aj — hereafter S93a; 
T93) that B — V colors provide reliable Tcff values at least 
for F and G-type stars, the latter ones being the objects 
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of main interest for our investigations; in fact, photomet- 
ric temperatures obtained from B — V color indices are 
in general in good agreement with temperatures obtained 
via spectroscopic analysis (see, e.g., Sestito et al. l2UU4f) . 
On the other hand, B — V are not very good T^g indi- 
cators for cool and/or active stars: first of all there is an 
intrinsic effect due to the fact that at lower Tes the peak 
of the blackbody distribution shifts towards lower ener- 
gies (higher wavelengths); for this reason, the dependence 
of the effective temperature on B — V becomes rather flat 
for B -V above - 1 (see, e.g., Grav lTFTSl Bessel lTOT^ . 
Second, as discussed by S93a, in ultra-fast rotators, which 
represent a significant fraction of young late-type objects, 
B — V colors can be altered either by effects of surface 
activity (e.g., spots; Stauffer et al. '198^ or the color dis- 
tortion may arise from structural changes due to rotation. 

In order to check on possible differences among Toff de- 
rived from different color indices, particularly for young 
cool and active stars, we show in Fig.^a comparison be- 
tween the Teff derived by us from BV photometry (see 
below) and those published in the original papers for the 
young cluster IC 2602 (Randich et al. [TM7I IMJTjl . derived 
as an average of the results from two calibrations (based 
on B — V and V ~ I). The figure indicates that the two 
sets of temperatures are in good agreement, not only for 
warm stars, but also for cooler objects, implying that the 
use oi B — V colors should not introduce major internal 
errors in the present study. 

Effective temperatures were computed by us from 
dereddened B — V colors using the calibration of S93a: 
Teff = 1808(B - y)o^-6103(B-T/)o-h8899 K; the authors 
found this temperature scale to be reasonably consistent 
with the temperatures used for the Pleiades by Boesgaard 
& Friel 1)199 01. who performed a calibration using F-type 
stars, and with that of Arribas & Rogers H1989|l derived 
from cool stars with the IR flux method. 

3.2. Li abundances 

The Li i feature can be blended with a nearby Fe i line 
at 6707.44 A, depending on the spectral resolution; for 
this reason, if necessary, the EWs were corrected from 
the Fe contribution following the prescriptions of S93b: 
EW{Fe)=20{B - V)o - 3 mA. For the following samples 
the resolving power was high enough to allow the separa- 
tion of the two features: Pleiades of Jones et al. (|1996|l . 
M 34 (Jones et al.T^, NGC 752 (Sestito et al. ^5007, 
Hobbs & Pilachowski J986a) and NGC 188 of Randich et 
al. I|2003(l . On the other hand, James & Jeffries l|1997(l and 
James et al. f^UOU' for NGC 6475), Jeffries & James iP^Hl 
Blanco 1) and Balachandran et al. ( 1988t 119961 a Per) 
published EWs already corrected from the Fe i contribu- 
tion but with a method different from that of S93b. Sestito 
et al. lPD03^) carried out a Li analysis of NGC 6475 and 
showed that, for stars in common with James & Jeffries 
l|1997|l and James et al. ()2000|l . there is a good agree- 
ment between the different sets of deblended EWs. Also 



Randich et al. H1998|l . in the analysis of a Per trusted 
in the Fe i correction performed by Balachandran et al. 
(flM^ [TMH| . Thus, we could safely use the deblended 
equivalent widths reported by these authors. We finally 
mention that upper limits in EW were never corrected 
for the Fe i blending. 

Li abundances were computed from T^s and deblended 
EWs by an interpolation of the curves of growth (COG) of 
S93b. These authors tested the accuracy of the results ob- 
tained from the COGs and reduced the possible systematic 
errors by calculating also synthetic spectra for broad-line 
stars. The two techniques were found to be in good agree- 
ment, while discrepant results for a subsample of stars 
were attributed to the saturation of the lines; a brief dis- 
cussion of systematic errors due to different methods of 
analysis is reported in Sect. 13. 3T^ 

The COGs of S93b are based on the assumption of 
local thermodynamic equilibrium (LTE), thus Li abun- 
dances were corrected for non-local thermodynamic equi- 
librium (iVLTE) effects — which are very important es- 
pecially in the case of cool stars — using the prescription 
of Carlsson et al. H1994|l . 

3.3. Uncertainties in Li abundances 
3.3.1. Random errors 

Random errors in Li abundances derive mainly from un- 
certainties in effective temperatures — due to photome- 
try uncertainties — and in EWs, while other parameters, 
such as surface gravities and microturbulence velocities do 
not significantly affect the determination of Li abundance 
from the Li i feature. 

Random errors in EWs come obviously from uncer- 
tainties in the measurement, depending on the S/N of 
the spectrum, and in addition, as mentioned in Sect. |21 
they depend on the spectral resolution, i.e. on whether 
or not the Li i doublet can be resolved from the nearby 
Fe I line. In order to check how large can be the error 
introduced in estimating EW{Fe) via the correction of 
S93b, we measured its value in spectra of stars with dif- 
ferent Toff and for which we were able to separate the 
two features. We found that for solar- type stars (Toff ~ 
5700 K; in NGC 752) the measured £;W^(Fe) range be- 
tween 9 and 12 mA, while one would obtain EW{Fe) 
~ 10 mA using the S93b formula; for hotter stars (Toff ~ 
6200 K) the observed EW{Fe) is 5-8 mA (analytical: ~ 7 
mA), while for cooler stars (in the Hyades) around 5000 
K the analytical estimate gives ~ 14 mA, to be compared 
with our measurement of EW{Fe) ranging from 13 up to 
18. Therefore, the largest difference between the measured 
EWs{Fe) and those estimated with the S93a prescriptions 
is of the order of 5 mA: this difference can be significant 
only in the worse case of rather cool stars (below 5000 
K) and low EW ('^ 10-20 mA) since it would correspond 
to an error in logn(Li) of ^ 0.15 dex. Note however that 
stars with Toff below 5000 K in young cluster have much 
larger Li EWs; on the other hand, the few measurements 
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at these temperatures in older clusters (e.g., Coma Ber) 
are upper limits in Li, from which the Fe i contribution is 
not subtracted. On the contrary, for solar-type stars in old 
clusters (e.g., with T^g ~ 6000 K and EW - 50 mA), even 
in the most conservative case, an error of 5 niA would cor- 
respond to 0.05 dex in logn(Li), which is well below the 
uncertainties usually estimated for these objects. 

Errors in EWs were published in the original paper 
only for a part of the clusters in the dataset; for these sam- 
ples we retrieved the published AEW and we assumed er- 
rors in effective temperatures according to the uncertain- 
ties in the photometry or in the same Teg, when quoted 
by the authors using Toff calibrations consistent with ours. 
Thus, we were able to estimate errors in Li abundances by 
quadratically adding AEW and ATcff- 

On the other hand, for about half of the sample (in- 
cluding the most populous ones — Pleiades of S93b and 
Jones et al. IT^ Hyades of T93 and Praesepe) AEW 
and ATeff for each star were not reported in the original 
papers and the authors quoted only typical average er- 
rors in these parameters or in log n (Li). In those cases we 
computed a mean error for each cluster starting from the 
typical average uncertainties. 

A detailed description of the error estimates for each 
cluster is given in Appendix. 

3.3.2. Systematic errors 

In order to investigate possible systematic effects on our 
temperature scale, in Figure |21 we compare the TcS val- 
ues determined by us for the 100 Myr old open cluster 
Blanco 1 with those derived by Jeffries & James fiyyyi): 
this is the sample for which we found the largest difference 
between our T^s scale and that used by previous authors. 
Jeffries & James ()1999|l derived effective temperatures for 
the hottest star {B — Vq < 0.64) from B — V colors using 
the calibration of Saxner & Hammarback H1985|l . which 
contains a metallicity dependent term; for cooler objects 
Jeffries & James applied the Bohm-Vitense (1981) so- 
lar metallicity calibration modified to take into account 
the higher than solar Fe content of this cluster ([Fe/H] 
= +0.14; Jeffries & James [T999|l . Random errors in our 
Toff were estimated from uncertainties in the photometry 
(ranging from A{B - T^) ~ 0.01 up to 0.03), while for the 
temperatures of Jeffries & James we considered an aver- 
age error bar similar to those plotted in their Fig. 1 (since 
they do not quote a AToff value in their tables). A reason- 
able agreement between the two samples is reached only 
for the two hottest stars; in the other cases our temper- 
atures are systematically lower, with differences of ~ 200 
K for stars around 5500 K, increasing up to almost 350 K 
for stars cooler than 5000 K. 

In the case of other clusters analyzed by the group 
of Jeffries (Coma Ber, NGC 2516, NGC 6633), all with 
slightly sub-solar Fe contents, the differences between our 
Toff and those quoted in the original papers, derived us- 
ing the metal-dependent Saxner & Hammarback (| 1985)1 



calibration, are rather small and within the errors: this 
indicates that the latter calibration is consistent with the 
S93a one, at least for [Fc/H]^0.0. The large discrepancies 
present in the case of Blanco 1 instead suggest that when 
over-solar [Fe/H] are considered, the metallicity correction 
for the Toff vs. B — V relationship might be important. 

As a test, we derived effective temperatures for 
Blanco 1 using also the metal dependent calibration by 
Alonso et al. H1999|l and assuming [Fe/H]=-|-0.14. We 
found that for stars with B — V < 0.70 the Toff derived 
with the Alonso calibration are cooler than ours by 50 
K. On the contrary, for larger colors, the T^s values found 
using the Alonso calibration are hotter than ours, with 
differences up to 250 K for the coolest star (Teff~ 4500 
K). This confirms that the [Fe/H] correction is not neg- 
ligible for over-solar metallicities and cool stars, and it 
represents an additional warning for the future users of 
the database. We stress however that only four clusters in 
the database have metallicities significantly higher than 
solar (Blanco 1, NGC 6475, the Hyades and IC 4651) and, 
as evidenced above, the Tcs-{B — V) relationship for solar- 
type stars — the objects of main interest for our discussion 
in Sect. 4 — is not strongly affected by [Fe/H]. 

Systematic errors for the Li distributions of open clus- 
ters are mainly due to the method of analysis; we de- 
rive Li abundances with the COGs of S93b, but slightly 
discrepant results would be obtained using a different 
method. In Table El we show a comparison between 
log7i(Li) values obtained with the COGs of S93b (Col. 5) 
and those obtained using the program MOOG (version 
2002 — SnedenHnZH Col. 4) for stars with various EWs 
and Toff in clusters of three different ages (NGC 6475, the 
Hyades and NGC 752). In the computation with MOOG 
we adopted the same surface gravity (logg — 4.5) for 
all the stars, while the microturbulence was calculated as 
g = 3 .2 X 10-''(Toff - 639 0)-1 .3(log g - 4.16) -h 1.7 ( Nissen 
nWU Boesgaard & Friel lTM)|l : we recall from Sect. lOH 
that Li abundances are not significantly affected by these 
two stellar parameters. 

Li abundances derived with MOOG and with the 
COGs appear to be in reasonable agreement: whereas val- 
ues derived with MOOG are systematically higher than 
the others, all the differences are < 0.1 dex, i.e. smaller 
than (or comparable with) typical errors quoted for the 
three clusters (see Appendix) ; note that the discrepancies 
in the final logn(Li) values decrease for lower Toff (in stars 
with similar EW). 

4. Results and discussion 

4.1. New Li abundances 

Figure |21 shows the comparison between the "new" Li 
abundances obtained with our method of analysis and 
those reported in the original papers for two open clus- 
ters: a Per (panel a) and Blanco 1 (panel b); filled sym- 
bols represent our estimates. Error bars represent indica- 
tive average errors in logn(Li) and Toff; in the case of a 
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Table 2. Comparison between Li abundances derived with the COGs of S93b and those derived with MOOG. 



Star 




EW 


log n(Li) 


log n(Li) 




(K) 


mk 


MOOG 


COGs (S93b) 


NGC 6475 










R14 


5888 


98 


2.85 


2.76 


R39A 


5696 


98 


2.64 


2.58 


R49A 


5772 


82 


2.61 


2.54 


R102 


5079 


67 


1.76 


1.69 


Hyades 










vB22 


5477 


56 


2.12 


2.03 


vB92 


5407 


15 


1.36 


1.34 


vB162 


5659 


59 


2.34 


2.25 


NGC 752 










P475 


5932 


59 


2.59 


2.50 


P520 


6151 


67 


2.86 


2.76 


P983 


5986 


45 


2.48 


2.40 



Per we found values very similar to those published in the 
original reference (for this reason we report only one error 
bar), while for Blanco 1 our errors are slightly more con- 
servative than those quoted by Jeffries & James 1)1999(1 . 

Data for a Per were taken from Randich et al. ([1998) ; 
note that we included in the plot also part of the sample 
of Balachandran et al. ifTS SS 1996) that was re-analyzed 
by Randich et al. As evident, differences in the Li distri- 
butions do exist, depending on the method of analysis. 
Figure OJi shows that in the case of a Per the two sets of 
Li abundances and Tes are characterized by random dif- 
ferences; however, in most cases the differences in Toff and 
logn(Li) are not very large and they lie within the errors. 

As shown in Fig. |3d and as already discussed, in the 
case of Blanco 1 there are significant systematic differences 
between our TcS and those retrieved from the original pa- 
per (see also Fig |2} ; this example shows how systematic 
errors due to different methods of analysis can affect the 
comparison between open clusters. Differences in Li abun- 
dances are only due to the discrepancies in Tcff since also 
Jeffries & James (|1999|l derived log n(Li) values using the 
COGs of S93b. However, the two sets of logn(Li) appear 
to be in good agreement, and the net effect is that the two 
Li distributions are shifted one with respect to the other. 

In summary Fig. O confirms that distinct methods of 
analysis result into random and systematic differences in 
the various sets of TcS and logn(Li); therefore the sign 
and the size of the error introduced by comparing clus- 
ters analyzed by different authors cannot be estimated a 
priori. This reinforces the conclusion that comparisons be- 
tween clusters of different ages and with stars covering a 
wide range of Tog greatly benefit from the use of a unique 
method of analysis. 

The complete dataset for Li in open clusters obtained 
with our method of analysis is presented in Figs.^Eland|ni 
in the first figure the Li distributions of young clusters 
(ages ranging from few Myr up to ~ 100 Myr) have been 
plotted, while in the other two we show the distributions 
for intermediate age (150-600 Myr) and old samples (^ 1 
Gyr.) Binary stars have been evidenced by filled symbols. 



We will not describe in detail the Li vs. Tcs distributions, 
whose features have been widely discussed in the original 
papers; we only note that the the star-to-star scatter is 
indeed evident both for cool stars in young clusters and 
among solar-type stars in the M 67, and that binary stars 
in general do not deviate from the average trends (see 
Sect.Ot. 

4.2. Time scales of Li depletion 

In Fig. 13 we plot the average Li abundance as a function 
of age for all the stars in the clusters in three tempera- 
ture ranges: Teg = 6200 ± 150 K (late F-type; panel a), 
Teff = 5900 ± 150 K (G-type; b) and Toff = 5600 ± 100 
K (late G-type; c); we chose a width of ±150 K for the 
first two intervals in agreement with the largest conserva- 
tive uncertainties in Toff for some clusters (see Appendix) ; 
in the case of stars around 5600 K we restricted instead 
the Toff range since in clusters as old (or older than) the 
Hyades the Li distributions are rather steep at these tem- 
peratures. Note that in some cases the cluster samples 
contain only few stars in the three Toff ranges considered; 
thus, in order to base our discussion on statistically more 
significant samples, the datasets of clusters with similar 
ages have been merged. Specifically, we merged the sam- 
ples of IC 2602 with those of IC 2391 and IC 4665 (age ~ 
30 Myr); a Per with NGC 2451 (~ 50 Myr); Blanco 1 
with the Pleiades (- 100 Myr); M 34 with M 35 and 
NGC 6475 (- 250 Myr); the Hyades with NGC 6633, 
Coma Ber and Praesepe (~ 600 Myr); the three clusters 
at 2 Gyr (NGC 752, NGC 3680 and IC 4651). Part of the 
merged samples include clusters with different Fe contents: 
however, as we have discussed in Sect. ^ (and as visible 
from Figs. 01 ElandEl), observations of cluster with the 
same age, but disparate [Fe/H] clearly indicate that overall 
metallicity does not affect Li depletion. For this reason, we 
can safely merge the datasets of clusters with similar age 
but different metallicity. We add again some words of cau- 
tion about this choice to remind that ages and [Fe/H] val- 
ues for the sample clusters are the results of independent. 
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and thus not homogeneous, investigations. Error bars in 
Fig. [7| indicate la standard deviations from the average; 
stars significantly deviating from the mean Li pattern in 
each merged sample have not been included; in the case 
of M 67 we computed the mean logn(Li) separately for 
the upper envelope and for the lower envelope. The av- 
erage values for the lower envelope in the 6200±150 K 
and 5900±150 K temperature ranges were computed us- 
ing ASURV that allows taking into account upper limits. 
For the coolest temperature interval we can instead only 
provide an upper limit to the mean, since we do not have 
enough data points. The average logri(Li) values for the 
various clusters (or groups of clusters) are summarized in 
Table El 

With only few exceptions, binaries have been included 
in the determination of the average abundance, since we 
find that their presence does not affect the Li distribu- 
tions, i.e. their Li abundances in general follow the av- 
erage trend for each cluster (Figs. ^ [S] and El see also 
Sestito et al. I21)U3|I : for example, among the few binaries 
excluded, there is star S1045 in M 67 (with T^g ^ 6100 K) 
since, as discussed in Pasquini et al. (1997,) . it lies signifi- 
cantly above the mean Li trend. Also star RX76A (Tog ^ 
4200 K) in NGC 2547 has an abnormally high Li abun- 
dance, and it is suspected to be a tidally locked binary 
in which Li depletion was inhibited (Jeffries et al. P003| ; 
anyway this is a very cool object, not included in our dis- 
cussion of the time scales of Li depletion. We mention in 
passing that Barrado y Navascues & Stauffer H1996|l in- 
vestigated a sample of Hyades binaries already observed 
by T93 finding that, whereas the general pattern is the 
same for double systems and single stars, binaries show 
slight over- abundances, that appear more conspicuous in 
tidally locked systems. Anyway, Barrado y Navascues & 
Stauffer obtained higher EWs for binaries with respect to 
T93; in our analysis we retain the EW values published by 
T93 and we find indeed that binary systems in the Hyades 
perfectly follow the average Li pattern of the cluster. 

The two dotted lines in the figure represent the allowed 
range for the initial Li abundance in Pop. i stars; the other 
curves represent models including non-standard physical 
processes (see below). Standard models are not shown in 
the plots since, as already mentioned, they do not predict 
any MS Li depletion; on the contrary, we are interested 
in the evolution of Li during this phase. As for the PMS 
phases, we recall that the amount of Li destroyed in this 
stage is a strong function of the assumed input physics 
in the models; different standard models result into dis- 
tinct amount of PMS depletion and thus into disparate 
logn(Li) values on the ZAMS (e.g., D'Antona & Mazzitelli 
HMM Siess et al. 2000"), but, as stressed in Sect.Q our aim 
is the investigation of the differential Li depletion during 
the MS, neglecting the absolute ZAMS abundances of the 
stars. 

Fig. Ell shows that late F-type stars (Tcff = 6200 ± 
150 K) deplete a certain amount of Li during the PMS 
phase; note that in this plot data for IC 2602, IC 2391 and 
IC 4665 are lacking (since no observations are available 



for stars within this Tee interval), thus a lower limit for 
Li abundances of ZAMS stars is represented by the mean 
logn(Li) of a Per -f NGC 2451 (age - 50 Myr). Then, 
only a small amount of destruction is suffered by these 
stars from 50 Myr to ~ 250 Myr; Li depletion slightly 
increases again beyond this age, as suggested from the 
points representative of clusters at the age of the Hyades 
and older (NGC 752 -I- IC 4651 + NGC 3680 and the 
upper envelope of M 67); Li destruction then appears to 
slow down between 2 and 5 Gyr (see below — plateau 
at old ages) . Finally, stars in the lower envelope of M 67 
have depleted an amount of Li larger by a factor ~ 4 with 
respect to their Li-rich counterpart. 

Panels b) and c) show that the Li evolution histories 
for stars in the Tes range 5500-6050 K are qualitatively 
similar to that of hotter stars: these objects deplete Li 
by a factor lower than 2 (depending on the assumed ini- 
tial content) during the PMS; then Li depletion virtually 
slows down/stops between the age of a Per up to ^ 150 
Myr in the case of panel c) and up to 250 Myr for panel 
b), and subsequently it accelerates again up to ~ 2 Gyr. 
After this age a plateau in average Li abundance is present 
in the case of the three oldest groups of clusters (i.e. up 
to the age of NGC 188), meaning that depletion might 
become ineffective, as discussed by Randich et al. H2003|l : 
note however that in panel c) only one star of NGC 188 is 
present, with logn(Li) much higher than that of other old 
clusters, thus it cannot be considered as significant in the 
overall comparison. In these T^s ranges the spread among 
stars in M 67 is even larger than for stars considered in 
panel a): the mean logn(Li) for the upper and lower en- 
velopes of M 67 differ indeed by a factor 8 for stars with 
Tcff — 5900 ± 150 K, and by a factor ~ 4 at least, for ob- 
jects with Tcff = 5600 ± 100 K. The Li content of the Sun 
is also shown in Fig.Eja and it is similar to that of the most 
Li-poor objects in M 67 (and in the solar neighborhood; 
Pasquini et al. ll994;i . 

Although stars in the three T^s ranges investigated are 
characterized by a similar qualitative behavior of Li deple- 
tion, it is important to note that the quantitative amount 
of surface Li destroyed within the various age intervals in- 
creases for cooler stars. The "final" abundances beyond 2 
Gyr have indeed values of ^ 2.5, 2.3 and 1.9 for stars in 
panel a), b), c) respectively (and also the average abun- 
dance of the lower envelope in M 67 decreases — 1.9, 1.3, 
and less). 

In summary, this discussion suggests that late F and 
G-type stars within the three ranges (Toff — 6200 ± 150 
K, Toff = 5900 ± 150 K and Toff = 5600 ± 100 K) have 
qualitatively similar Li depletion histories which can be 
schemed as follows: 

— a small amount of PMS Li destruction; 

— plateau during the early stages of the MS; 

— MS depletion; 

— plateau at old ages. 

However, the cooler is the temperature, the higher is 
the amount of Li destruction suffered during the vari- 
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Table 3. Average logn(Li) in three ranges of Tea. 



dustcr(s) 


Age 


log 7i(Li) icT 


log 7i(Li)i(7 


log 7i(Li)iCT 




(Gyr) 












conn _L 1 c: n Tj^ 

-^off — 6200 ± lou K 


Jeff — 5900 ± 150 K 


Jeff — 5600 ± 100 K 


NGC 2264 


0.005 


3.34±0.10 


3.21±0.17 


3.18 


IC 2391/IC 2602/IC 4665 


0.03 




3.06±0.13 


2.90±0.11 


a Pcr/NGC 2451 


0.05 


2.97±0.13 


2.96±0.10 


2.94±0.09 


Plciadcs/Blanco 1 


0.1 


2.95±0.09 


2.96±0.07 


2.79±0.10 


NGC 2516 


0.15 


2.92 


2.85±0.11 


2.80±0.23 


M 34/M 35/NGC 6475 


0.25 


2.92±0.13 


2.79±0.14 


2.57±0.19 


Hyades/Pracscpe/Coma Bci /NGC 6633 


0.6 


2.77±0.21 


2.58±0.15 


2.14±0.22 


NGC 752/NGC 3680/IC 4651 


2 


2.65±0.13 


2.33±0.17 


1.91±0.24 


M 67 — upper envelope 


5 


2.55±0.18 


2.25±0.12 


1.90±0.15 


M 67 — lower envelope 


5 


1.92±0.05 


1.32±0.09 


<1.29 


NGC 188 


8 




2.34±0.14 


2.30 



ous phases and the shorter is the epoch of the early-MS 
plateau. 

Two features are worth being stressed: first, Li deple- 
tion is not a continuous process that can described with 
a t~" law. Indeed Li depletion is characterized by dif- 
ferent timescales in the various age ranges. Second, the 
current homogeneous analysis of Li abundances in open 
clusters covering a wide range of ages allows us to con- 
strain for the first time the age at which extra-mixing 
processes or other non-standard mechanisms start work- 
ing in solar-type stars and objects of similar temperatures. 
Whereas solar-type stars do not undergo further Li deple- 
tion during the early stages of the MS (convection is not 
sufiicient, since the convective zones — CZ — are too shal- 
low at these ages), an additional mechanisms driving Li 
depletion appears around an age of 150-250 Myr. 

4.3. Comparison with the models 

Several non-standard processes have been proposed in 
the last years to be included in models: mixing driven 
by rotation or by internal waves (see below); magnetic 
fields, which might either enforce rigid rotation in stars 
or produce differential rotation (e.g., Charbonneau & 
McGregor [rnnmnSllinS^; microscopic diffusion, deter- 
mining the separation of atomic species in the stellar gas 
(e.g.. Richer & Michaud [T?M1 Chaboyer et al. I1995ajl : 
mass loss through stellar winds (e.g., Swenson & Faulkner 
, causing Li dilution in stars; tachocline mixing (Piau 
et al. l2UU3fl . The tachocline is a thin layer of separa- 
tion between the CZ, in which differential rotation occurs, 
and the internal radiative core, rotating as a rigid body; 
tachocline mixing is a hydrodynamic process, not related 
to magnetic effects. 

We focus here on those mechanisms for which detailed 
predictions of Li depletion for solar-type stars as a func- 
tion of age and T^s are available. Namely: (i) slow mixing 
induced by rotational instabilities and angular momen- 
tum loss (AML; Chaboyer et al. ll905al[T995bl — hereafter 
CDP95b, dashed line in Fig.[7| Deliyannis & Pinsonneault 
119971 — DP97, dot-dashed line/cross in panel c)) and (ii) 
mixing driven by internal waves (Montalban & Schatzman 
IMini — MSOO, sohd fine). Note that the latter mod- 
els include only gravitational waves, neglecting rotational 



mixing. It is worth mentioning that very recently other 
models have been developed which take into account the 
transport of angular momentum through internal gravita- 
tional waves together with the effect of rotation (Talon & 
CharbonnelEOnSI)- 

In the models by MSOO (see also references therein) 
chemical mixing is induced by the non-adiabatic propa- 
gation in the stellar interiors of internal waves generated 
at the boundary of the CZ. Panel a) shows that these 
models do not predict any surface Li destruction in stars 
hotter than 6000 K during the whole age interval consid- 
ered, at variance with observations; on the other hand, too 
high Li abundances are predicted for cooler stars (panel 

b) ) up to an age of 3 Gyr. After this epoch the deple- 
tion suddenly increases and the theoretical value for stars 
with solar age is only slightly lower than that observed for 
the upper envelope in M 67; however, although models at 
older ages are not present, it is clear that the plateau for 
objects older than 2 Gyr is not reproduced and that the 
theoretical behavior is opposite to the empirical one. The 
same happens for objects with Tcs= 5600 ± 100 K (panel 

c) , but in this case the depletion predicted by the models 
strongly increases around 1 Gyr, due to the deeper con- 
vective zones of these stars. In conclusion, MSOO models 
are not able to reproduce quantitatively and qualitatively 
the observed time scales of Li depletion, thus mixing in- 
duced only by gravity waves can be excluded as the main 
process responsible for MS Li destruction. 

Rotation can drive Li depletion since meridional cir- 
culation and/or instabilities triggered by differential ro- 
tation in the stellar structure induce a slow mixing over 
time scales comparable to those of light element burning. 
The agent responsible for the rising of differential rotation 
would be an external torque, i.e. AML (spin down): this 
process acts to generate radial gradients of angular ve- 
locity and instabilities leading to AM redistribution and 
mixing of material. Stars born as fast rotators can either 
undergo a large amount of AML during the PMS or the 
early stages of the MS and lose a large amount of their 
initial Li content, or they can continue to rotate rapidly 
preserving a high log n (Li); initial slow rotators cannot 
suffer a large amount of AML and as a consequence the 
stellar structure adjusts itself on a circulation-free state. 
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We report in the three panels of Fig. [7| two different 
rotational models. The models by DP97 include AML and 
AM transport dominated by the secular shear instability 
fZahn ni)92|l . while other parameters such as the velocity 
estimates for meridional circulation and other instabilities 
were not varied in that study (see also Pinsonneault et 
al. ll989[ 119901 for a detailed description of the models). 

The CDP95b models (see the quoted reference and 
Chaboyer et al. I1995a|l include both AML and meridional 
circulation, treated as a diffusive process. In addition these 
authors take into account also diffusive mechanisms not 
related to rotation (e.g., gravitational settling). Chaboyer 
and collaborators computed several models by varying a 
series of parameters, namely: the mixing length parame- 
ter; the He and metal contents {Y and Z); the efficiency 
of overshooting; a constant multiplying the diffusion coef- 
ficients (/m); and five rotational parameters including the 
initial rotational velocity and a critical angular velocity at 
which AML reaches a saturation level. We chose to repre- 
sent the "best fit" model by CDP95b, i.e. the one which 
in their opinion provides the best agreement with obser- 
vations of open clusters. This model is labeled as VN in 
the reference paper; the initial rotational velocity is 30 km 
s-i. 

Figure shows that in some cases it is present a rea- 
sonable agreement between the predictions of rotational 
models and the empirical patterns. To be noticed that the 
models by CDP95b indeed predict a discontinuous Li de- 
pletion. More specifically, for stars with Tee= 6200±150 K 
(panel a), the model by CDP95b gives a rather good qual- 
itative agreement with the empirical behavior — at least 
up to 2 Gyr; whereas the theoretical early MS plateau is 
too high with respect to the empirical one (no PMS deple- 
tion is predicted), as already stressed, we are only inter- 
ested in the differential amount of MS Li depletion, thus 
the absolute Li value in the plateau is not important. On 
the other hand, the slowing down of Li destruction after 
2 Gyr is not reproduced adequately. The DP97 isochrones 
are apparently not able to well fit the shape of the early 
MS plateau, but actually models are available for three 
ages only (100 Myr, 1.7 Gyr and 4 Gyr). Thus, data be- 
tween the early MS and the age of the Hyades are lacking 
for these models and we cannot draw any conclusion. In 
any case, the DP97 model with Wrot — 30 km s^^ is in 
strident quantitative disagreement with the empirical ev- 
idence. 

Similar arguments apply to stars in panel b) (Toff= 
5900 ± 150 K) and c) {Tcs= 5600 ± 100 K). The two pan- 
els show that the value of the early MS plateau is only 
slightly overestimated by the GDP95b models, while a 
good agreement is present between the Hyades age and 
2 Gyr. However, the plateau for old stars is not fitted, i.e. 
Li abundances decrease towards the solar age. Also the 
models by DP97 are not able to reproduce the shape of 
the late MS plateau. A too strong depletion is predicted by 
the model with Vi-ot = 30 km s^^ and Tcff = 5900 K, if com- 
pared with open clusters; this track would be able instead 
to fit the solar abundance; only one logn(Li) value (at the 



Pleiades age) is available for the track with Vrot — 30 km 
s~^ and Tctf— 5600 K (represented by a cross in panel c)), 
much lower than the empirical one. 

In summary, the approach of models including slow 
mixing related to rotation (and diffusion) appears reason- 
able since these models can somehow reproduce the ob- 
served Li depletion pattern. The CDP95b tracks are able 
to reproduce the discontinuous rate of MS Li depletion up 
to the Hyades age, while no conclusion can be drawn for 
the DP97 models. 

The models by CDP95b are probably best suited for 
reproducing the early-MS plateau and the subsequent ac- 
celeration of Li depletion, due to the competition between 
rotation and diffusion. The theories by DP97 ignore mi- 
croscopic diffusion, which instead can have an important 
role in particle transport: for example (see Chaboyer et 
al. E95a) the diffusion of ^He induces gradients in the 
mean molecular weight, which inhibit rotational mixing 
(since extra energy would be needed to move material to 
zones of different molecular weight). On the opposite, ro- 
tational mixing might inhibit the separation processes of 
diffusion. In this sense, the plateau after the ZAMS could 
be the result of this competition. Also AML and merid- 
ional circulation might be competitive mechanisms; such 
combined effects are not present in DP97 models. As fi- 
nal comments, we remark that neither the CDP95b mod- 
els nor those of DP97 can reproduce the plateau at old 
ages. In addition, these rotational models are not able to 
reproduce the solar rotational profile measured by helio- 
seismology. In a very recent paper. Talon & Charbonnel 
|2005) presented models including gravitational waves and 
rotation, showing that the combined effect of such mech- 
anisms leads to Li abundances for stars at the Hyades age 
in agreement with observations, and at the same time they 
are able to shape the solar rotational profile. 

5. Conclusions 

We presented a new homogeneous analysis of all the Li 
data available in the literature for F, G and K-type stars 
in open clusters. The need of a dataset with all the clus- 
ters on the same scale of Li abundance and Tcff arises 
from the fact that comparisons between samples analyzed 
with different methods are affected by random and sys- 
tematic errors which cannot be estimated a priori. We 
have indeed shown that differences between the methods 
of analysis do exist. The database will be made available 
on the web, in order to allow homogeneous comparisons 
between the various samples and to test updated models: 
it includes 21 open clusters observed by various authors 
spanning the age range 5 Myr - 8 Gyr and with metallic- 
ities -0.2;^[Fe/H] ^-1-0.2. 

Based on the new Li abundances obtained by us for 
open clusters, we have determined average Li abundances 
for stars in different temperature bins and age ranges in 
order to investigate the time scales of Li depletion from 
ZAMS to late MS. For the first time we were able to put 
a limit in age to the appearance of non-standard mix- 
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ing mechanisms at work on the MS: stars do not appear 
to destroy Li in the early stages of MS, but Li depletion 
starts again around an age of ~ 150-250 Myr. Then the 
destruction is ineffective beyond ~ 1-2 Gyr in late F and 
G stars, except for a fraction of stars in M 67. Overall 
we showed that MS Li depletion is not a continuous pro- 
cess. We also confirmed the existence of a plateau at old 
ages. A caveat to our results is represented by the fact 
that ages and metallicities for the clusters have not been 
derived with a homogeneous method. Whereas our quan- 
titative estimates of time scales may be somewhat affected 
by this fact, we believe that our qualitative results, and 
in particular, the existence of a Li plateau at old ages do 
not change. 

As to the possible extra-mixing mechanisms operating 
during the MS, we compared the empirical scenario with 
the predictions of models including slow mixing related 
to gravity waves (Montalban & Schatzmann 2000') and 
to rotation and angular momentum loss (Chaboyer et al. 
I1995bl Deliyannis & Pinsonneault [TM7|l . 

We found that gravity waves (at least when not cou- 
pled with rotation — models by MSOO) can be excluded 
as the main agent responsible for MS Li depletion, while 
slow mixing induced by rotation might explain to some ex- 
tent the empirical behavior. In particular, the models by 
CDP95b including AML, meridional circulation and dif- 
fusion are those which appears at the moment to be the 
most appropriate. Nevertheless, several uncertainties are 
still present and none of the models proposed is able to 
reproduce the Li plateau observed for solar-type stars in 
clusters older than the Hyades. Understanding this empir- 
ical behavior is instead very important for its connection 
with Primordial Nucleosynthesis and Galactic evolution. 
We conclude that improvements in theoretical work are 
certainly needed for a proper comprehension of the prop- 
erties of mixing and Li evolution in stars. 
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Appendix 

We report here a more detailed description of the com- 
putation of random errors in Li abundances (due to un- 
certainties in Tqs and EW) for the open clusters in- 
cluded in the database and presented in Table ^ (see also 
Sect.ESUJ). 

The clusters with published errors in EWs are the fol- 
lowing: 

IC 2602: Randich et al. ipUTjl published the AT^b for 
each star, but their temperatures were derived in most 
cases from both B — V (with the S93a calibration) and 
V ~ I colors, thus in principle we cannot use directly 
these values. Since no error is quoted in the UBV pho- 
tometry used (Prosser et al. I1996|l , we adopted a con- 



servative ATcff = ±100 K, as Randich et al. H2001|l did 
for stars with Tcs derived from only one color index. 

IC 2391: the same as for IC 2602 (Randich et aL lM7T|l . 

IC 2547: Jeffries et al. ij^UUS^ give errors in T^g of ± 
100 K, using V - I colors from Naylor et al. H2002fl . 
Errors in B — V from the same photometric source can 
be as large as 0.03, corresponding to the same value 
ATeff - ±100 K. 

a Per: for the sample of Randich et al. H1998|l we used 
the ATeff quoted in the paper (~ ± 100, 150 and 200 
K), since no error information is present in the source 
of photometry (Prosser & Randich 1998 Prosser et 
al. riSgSfl : Randich et al. derived T^s by averaging two 
different calibrations (Bessel 119791 Alonso et al. I1996|l 
which provide temperature errors very similar to those 
obtained with the S93a calibration. For Balachandran 
et al. ifTOSSI ITMSjl see below. 

NGC 2451: we used the ATes of ±200 K quoted by 
Hiinsch et al. I|2004|l . which derived the temperatures 
from B — V colors with the S93a calibration. 

Pleiades: (i) Garcia Lopez et al. H1994|l : the same as for 
a Per of Randich et al. 119981 (but Garcia Lopez et 
al. use three different calibrations); (ii) Jeffries H1999b|l 
does not compute neither Li abundances, nor the tem- 
peratures (he only measures the EWs); we adopted 
ATeff — ±100 K, which is a conservative value for 
Pleiades stars. 

Blanco 1: we derived ATeff from errors in the photometry 
(0.01 <AB-V< 0.03). 

NGC 2516: we derived ATeff 70 K from errors in the 
photometry. 

M 35: we derived ATeff ^ 50 K from errors in the pho- 
tometry. 

NGC 6475: Alogri(Li) were already derived by Sestito 
et al. ()2()()3|l . also for the samples of James & Jeffries 
H1997|l and James et al. H2[)()()|l . using typical conser- 
vative ATeff of ± 100 K. 

Hyades: we were able to compute the errors for the two 
stars of Soderblom et al. (1990) for which we assumed 
ATeff~ 50 K (A(B - V) ~ 0.01). 

Coma Ber: in the three source papers the authors quoted 
the uncertainties in the photometry; we used these val- 
ues to compute the errors in Teff (ranging from 25 to 
150 K). 

NGC 6633: conservative ATcs derived from errors in the 
photometry are of ~ ± 100 K. 

NGC 752: Sestito et al. pTOjl computed errors for their 
sample and for that of Hobbs & Pilachowski fl9 86a|l 
assuming errors in Teff of ± 150 K (deriving from 
A{B — V), see the discussion in the quoted reference). 

IC 4651: ATcff = ±50 K from the photometry, with the 
exception of two stars which have AT^s = ±100 K. 

NGC 3680: we computed Alogn(Li) for the two stars in 
the sample of Randich et al. H2000|l using their AEW 
and assuming ATog = ±100 K as in Pasquini et al. 
1)20011 see below — errors in colors are not quoted in 
the photometry source by Nordstrom et al. I1997|l . 
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M 67: only Randich et al. H2Ul)2|l give errors in EWs. For 
their data, we considered AT^e = ±50 K, as in Jones 
et al. l(TM^ see below). 

NGC 188: Randich et al. 1)2003(1 assumed conservative 
uncertainties in T^g of ± 100 K deriving from the pho- 
tometry. 

These are instead the clusters for which we roughly 
estimated an average uncertainty in Li abundances, since 
errors in EW are not quoted: 

NGC 2264: (i) Soderblom et al. (TM^ quote only a mean 
error in Li abundance, Alogn(Li) = 0.1 dex; we re- 
tained this value, since we used the same method of 
analysis, (ii) King H1998|l give uncertainties of about 
20-30 mA in EW and of ±140 K in T^s for his sam- 
ple. From these values he obtained Alogn(Li) ~ 0.25, 
in agreement with us. Note that King used a different 
Teff calibration fBessel lT979|l however perfectly consis- 
tent with that of S93a. 

IC 4665 (Martm & Monies ITMTjl : AEW^IB mA, 
ATcff = ±150 K, using the calibration of Arribas & 
Martinez Roger H1988|l : we found this calibration to 
give errors in effective temperature very similar to 
ours. Martin & Monies quote an average final error in 
Li Alogn(Li) 0.3, but we find slightly lower values, 
Alogn(Li) - 0.23. 

a Per: Balachandran et al. 1)19881 11996)1 quote 
A log ri(Li) = ±0.1 deriving from uncertainties in EWs 
and Alogn(Li) = ±0.2 from uncertainties in Toff of 
±200 K. We cannot use this latter value, since, besides 
the different calibration, we adopted a different pho- 
tometric source (Prosser & Randich 119981 Prosser et 
al. 11998)1 ■ where errors in i? — are not given; as to the 
error deriving from EW, we can trust the value quoted 
by Balachandran et al. since they derived Li abun- 
dances with MOOG, which appears to be consistent 
with the COGs of S93b (see Sect.EXH. Thus, we can 
compute indicative errors by assuming ATeff ~ ±100 
K (quoted for most stars by Randich et al. 11998)1 and 
Alogri(Li) = ±0.1 deriving from EW uncertainties: 
we obtain errors in logn(Li) ranging from ^ 0.13 for 
the warmest stars up to ^ 0.17 for the latest spectral- 
type objects. 

Pleiades: (i) S93b give AEW—15 mA, without estimat- 
ing the uncertainty in Tog. Considering AT^s — ±100 
K, which we think to be conservative for this well stud- 
ied cluster, we find average Alogn(Li) ^ 0.13 — 0.17 
(the error increases for the coolest stars), (ii) Jones 
et al. 1)1996)1 : A£'W^=5 mA for slow rotators and 
AEW=20 mA for fast rotators. Assuming ATcff = 
±100 K we obtain a conservative Alogn(Li) of 0.15 
dex, both for fast rotators (which are mainly Li rich 
stars) and for slow rotators, (iii) Butler et al. 1)1987)) : 
they quoted an uncertainty of ^ 15 % in EW; as for 
the Teff, we used the errors quoted by Garcia Lopez et 
al. 1)1994)1 . who reanalyzed the Butler sample; errors in 
logn(Li) range from 0.15 to 0.40 dex. 



M 34 (Jones et aL [TM7|l : ATeg = ±100 - 130 K (using 
the calibration of S93a) and relative EW uncertainties 
of about 10%. It follows Alogn(Li) - 0.13 - 0.15. 

Hyades: (i) T93: AEW=2 mA, ATeff ±50 K, 
from uncertainties in B — V of ^ 0.01 (conserva- 
tive; see the detailed analysis in the original paper); 
we find Alogn(Li) = 0.1 for cool Li poor stars and 
A log n(Li) = 0.05 for hotter stars. 

Praesepe (Soderblom et al. I1993c)l : AEW^5 mA, no 
Teff uncertainty quoted. Assuming ATeff — ±100 K, 
errors in logn(Li) range from 0.1 dex for late F and 
G-type stars up to 0.17 dex for later spectral types. 

NGC 3680: for stars of Pasquini et al. ij^OTTIl . we started 
from ATcff — ±100 K quoted by them and we assumed 
a conservative uncertainty of 15% in EW. Errors in 
logn(Li) turned out to be ~ 0.10-0.15 dex. Note that 
they used the calibration of Alonso et al. 1)1996)1 , which 
results into errors similar to ATcff deriving from S93a. 

M 67: Jones et al. lPM?|l quote average AEW of ±10 
mA and ATeff = ±50 K for his data (using the S93a 
calibration), as well as for the previous samples (see 
Table ^1 re-analyzed by them. These values lead to 
Alogn(Li) 0.06 — 0.08 for stars with Li detection. 
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Fig. 1. Comparison between effective temperatures derived by us using the calibration of S93a and those presented in 
the original papers for IC 2602 (Randich et al. .1997. ■2001;) . computed by averaging two cahbrations based on B — V 
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Fig. 2. Comparison between effective temperatures derived by us using the calibration of S93a and those presented 
in the original paper for Blanco 1 (Jeffries & James IT999|I . computed with the calibrations of Saxner & Hammarback 
l(TU5^ and of Bohm-Vitense ipMT|l . 
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Fig. 5. Li distributions for intermediate age clusters (~150-600 Myr). 
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Fig. 6. Li distributions for old open clusters (age ^1 Gyr). 
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